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Observational constraints on loop quantum cosmology 
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In the inflationary scenario of loop quantum cosmology in the presence of inverse- volume correc- 
tions, we give analytic formulas for the power spectra of scalar and tensor perturbations convenient 
to compare with observations. Since inverse-volume corrections can provide strong contributions to 
the running spectral indices, inclusion of terms higher than the second-order runnings in the power 
spectra is crucially important. Using the recent data of cosmic microwave background and other 
cosmological experiments, we place bounds on the quantum corrections. 
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One of the motivations to search for a quantum the- 
ory of gravity is the desire to unify general relativity 
with quantum mechanics and, thereby, resolve classical 
singularities such as the big bang or those associated 
with black holes. Observational implications of quan- 
tum gravity, however, present a delicate issue. Based 
on dimensional grounds, cosmology in a nearly isotropic 
setting seems to allow quantum corrections only as pow- 
ers of the small quantity i-p\H « 10^^", where ip\ is the 
Planck length and = a /a is the Hubble radius (a is 
the scale factor in the fiat Friedmann-Robertson- Walker 
background and dots denote derivatives with respect to 
cosmic time t). This dimensional argument is supported 
by low-energy effective actions of higher-curvature type. 

Dimensional arguments, generally, are overcome if 
there are more than two dynamical scales of the same 
dimension. Detailed physics rather than rough estimates 
are then required to determine which geometric mean of 
the scales is relevant in a given regime. In cosmology, an 
additional distance scale L would allow a multitude of 
dimensionless combinations ip^H^L'' with a — /3 + 7 = 0, 
not all of them small. Quantum gravity provides ample 
motivation for the existence of a third scale by suggesting 
discrete spatial structures. While the discreteness scale 
L is often expected to be near £pi, it is not identical to 
it and also depends on excitation levels of states (rather 
than just Newton's and Planck's constants). 

One explicit formulation of such a discrete version of 
gravity is loop quantum gravity (LQG) [ij. Discrete- 
ness arises on the space of metrics (geometrical oper- 
ators acquiring discrete spectra). In a nearly homoge- 
neous quantum space-time, one can think of any region 
of volume V to consist of discrete patches, each roughly 
of size with the length L determined by an underly- 
ing quantum-gravity state. Discrete spectra imply that 
derivatives by L, as they ubiquitously appear in canon- 
ical expressions via Poisson brackets, are replaced by fi- 
nite difference quotients. As a simple example for so- 



called inverse-volume corrections, {2\/L)^^ ~ d\fL/dL 
would, when evaluated for discrete operators, become 
{\/L + ^pi — \/L ~ £pi)/2£pi, which strongly differs from 
{2\/L)~'^ for L ^ ^pi. For larger L, corrections are per- 
turbative and of the order £pi/L; no factor of H appears. 
The ratio £pi/L can easily be much larger than ipiH, ex- 
plaining why this type of discreteness could give rise to 
stronger quantum effects. 

The results of detailed constructions in LQG, following 
0: 01' ^^11 be summarized momentarily. First, we em- 
phasize that the discreteness does not break general co- 
variance in the equations used here (assuming small cor- 
rections). This has been demonstrated by an elaborate 
analysis of the gauge contents of the quantum-corrected 
theory, verifying the existence of a closed algebra of gauge 
generators Covariance, and the space-time structure 
it belongs to, is then not destroyed but deformed. (Defor- 
mations of classical symmetries play an important role in 
several approaches to quantum-gravity phenomenology 
0. The deformations considered here are on a differ- 
ent footing, however, because they do not refer just to 
Poincare transformations of Minkowski space.) 

Here, using currently available data, we place con- 
straints on inverse- volume corrections for inflation. Since 
scalar and tensor perturbations are subject to strong 
modifications of the power on large scales, the corrections 
are bounded from above. A detection of gravitational 
waves and the precise measurement of cosmic microwave 
background (CMB) anisotropics in future observations 
such as Planck will potentially allow us to make a deci- 
sive test for loop quantum cosmology (LQC) inflation. 

A simplified implementation of corrections expected 
from LQG in cosmological scenarios via perturbations 
around homogeneous or other reduced models can be 
achieved in LQC Q . With a phenomenological approach 
to effective dynamics, the cosmological equations can be 
summarized in a single Mukhanov equation for the gauge- 
invariant scalar perturbation Uk, u'^-^- {s'^ k'^ ~ z" / z)uk = 



2 



in momentum space with the comoving wave num- 
ber fc, where primes denote derivatives with respect to 
conformal time t = J dt. Similarly, tensor modes 
are subject to the equation w'^ + {a^k'^ — a" /a)wk = 
Here, z(a, i^) and a{a) are background functions and 
« 1 -I- 2ao(5pi and = 1 + xf^pi arc the propaga- 
tion speeds squared, differing from the speed of light by 
quantum corrections. 

The quantum corrections are characterized by (i) nu- 
merical coefficients ao and x and (ii) the function Jpi (x 
a~'^ determining the size of inverse-volume corrections. 
The values of ao, x and a are currently subject to 
quantization ambiguities, x is parametrized as x = 
(jvq[(j/& -f 1)/3 -I- ao(5 — cr/3)/2, where is related to 
ao and a by the consistency condition 



Vq{(j — i) ~ 3ao (cr — 6)/ (cr -I- 6) . 



(1) 



While tJ takes values in the range < fi < 6, the size of 
(5pi docs not depend on the values of ao and vq. With 
cr > 0, (5pi is larger at early times, in agreement with 
discreteness departing from the Planck scale in a more 
classical universe. The aim of this Letter is to restrict 
(5pi by observations. We will mainly place bounds on 
the combination ao^pi during slow- roll (SR) inflation, for 
which the precise origin of ao and v^) or the scale hidden 
in 5pi is not essential. 

Corrections in the evolution equations arise only in the 
fc^ term, not in the time derivative of the d'Alembertian, 
yet they are covariant according to the corrected gauge 
transformations Thus, one typical assumption of 

higher-curvature theories is violated. Moreover, the 
propagation speed of tensor modes differs from the scalar 
one since in general 2ao ^ X- Again, this is only possible 
with the change in the underlying manifold and gauge 
structure, and gives rise to additional characteristic ef- 
fects. With different types of equations for scalar and 
tensor modes, there are changes to the standard infla- 
tionary spectra and the tensor-to-scalar ratio. 

In Ref . [1] , two of us evaluated the inflationary observ- 
ables in terms of the three SR parameters e = —H/H^, 
T] = —(f/{Hip), and = (if/ip)/!!^, where is a scalar 
field with potential V{(p). In order to place observa- 
tional bounds on concrete inflaton potentials, it is more 
convenient to use SR parameters expressed by V and 



its derivatives: 

I = K~^V,ipV,^^^/V^ where 
tational constant 
to e. 



SttG (G is the gravi- 
For conversion formulas from e, ?7^f ^ 



'iVviiv and all the technical details we refer to 



together with a discussion of cosmic variance. 

The power spectra of scalar and tensor perturbations, 
evaluated at the Hubble horizon crossing during inflation 
{k w clH), are given, respectively, by Q 



7re 



(l + 7td>i), (2) 



where 7s = i'o(cr/6-|-l)-|-(Tao/(2e)— [^7l/o(c^-|-6)-|-3ao(15 — 
cr)]/[18(cr -h 1)] and 7t = (cr - l)ao/(cr + 1). We expand 
the scalar spectrum about a pivot wave number fco, as 



Inn(fc) = l^Vs{ko) + [ns{ko) -l]x 



as{ko) 



m— 3 



(3) 



where x ~ ln(fc/fco), ns{k) — 1 = dln'Ps{k)/dlnk, and 
ai™\k) = d™~^as/(dlnfc)™^^. The tensor spectrum 
can be expanded in a similar way with a different in- 
dex nt(fc) = d\n'Pt{k)/d\nk. While such expansions to 
second order are standard in cosmology, terms of order 
higher than 2 will become important in our analysis. 
The spectral indices are 



1 



PI : 



-2ev-c„,(5pi, (4) 



with quantum-gravity corrections ^ = /s,t + • ■ ■ whose 
dominant contributions are fs = cr[3ao(13cr — 3) + t'ocr(6 + 
llcr)]/[18(cr + 1)] and ft = 2a'^ao/{(T + 1). For a > 0(1) 
the variation of i5pi is fast {Spi cx a^'^ cx k^"' at Hub- 
ble crossing), so that fs^t provide dominant contribu- 
tions to the scalar and tensor runnings as well, as,t(fco) = 
dns^t/dlnklk^ka ~ o'/s,t'5pi(fco)- Similarly, the mth order 
terms are al"^\ko) sa (— l)™cr™~-'^/s t(5pi(fco) and hence 
we can evaluate the sum in Eq. ^ as 
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fs,tSp\. (5) 



This expression is valid for any value of cr and of the pivot 
scale fco within the observational range of CMB. Since 
the LQC corrections to the runnings ag.t can be large, 
inclusion of the higher-order terms ^ is important to 
estimate the power spectra properly. 

For the CMB likelihood analysis we also take into ac- 
count the second-order terms of slow-roll parameters, 
i.e. as — — 24e^ -I- IGevrjv — 2^^ -I- Ca^Spi and at = 
— 4ev (2ei, — 77v) -l- Cq,j (Spi, where the dominant contri- 
butions to Cq,^ J correspond to Cq,^ ^ « cr/s,t. In the 
numerical code, the full expressions of the coefficients 
c„^ ^ and ^ [8|] are used. At the pivot scale ko we 
have the tensor-to-scalar ratio r(fco) = 'Pt{ko)/'Ps{ko) = 
16ev,(fco) -I- CrSp\{ko), where c,, = 8[3ao(3 -I- 5cr -I- 6cr^) — 
z/o't(6 + ll(T)]ev(A:o)/[9((T + 1)] - 16(Tao77v(fco)/3. 

In the quasi-de Sitter background, Spi cx k~"' gives 
the relation Spi{k) « Spi{ko){k/ko)-'' = Spi{io){i/io)~'' , 
where £ are the CMB multipoles related to k via k « 
{h/10^)£ Mpc-i {h « 0.7 is the reduced Hubble con- 
stant). With the large- volume expansion of quantum 
corrections, we require that Sp\{k) ^ 1 at all scales. 
For cr > the LQC correction is most significant on the 
largest scales observed in the CMB {£ = 2). This prop- 
erty can be clearly seen in Fig. [U where the pivot scale 
for the scalar power spectrum is taken to be £o = 29. In- 
tuitively, this happens because the largest cosmological 
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FIG. 1. Primordial scalar power spectrum 'Ps(^) for the case 
n = 2, (T = 2, and ev(fco) = 0.009 with three different values 
of (5(fco): (classical case, dotted line), 7x 10~^ (experimental 
upper bound, solid line), 4.8 x 10"* (1/10 of the a-priori upper 
bound, dashed line). Here the pivot wave number is ko — 
0.002 Mpc~\ which corresponds to io = 29. 



scales correspond to those where and when space-time 
quantum effects were larger, while smaller scales have 
been affected by ordinary physics. Imposing the condi- 
tion 6pi{£ = 2) <C 1, this gives 



5pi(4) « {2/eor 



(6) 



at the multipole £o- For larger a and £o, 6p\{io) is con- 
strained to be smaller. We assume that inflation com- 
pletely takes place in the large-volume regime, and that 
longer-wavelength modes that might violate the bound 
^ do not enter the inflationary Fourier analysis. 

For concreteness, let us consider the power-law poten- 
tial V{ip) = A(p", for which = n^/(2K^(/3^) and 



2(n-l) ,2 _ 4(n-l)(n-2) , 



Among the variables a, ao, and vq we have the rela- 
tion ([1]), a condition under which, for given n and cr, 
the inflationary observables can be expressed via and 
S = aoSpi for cr 7^ 3, or by e,,- and 6 = I'o^pi for cr = 3. 

We carry out the CMB likelihood analysis by varying 
the parameters ev and 6 in the cosmological Monte Carlo 
(CosmoMC) code ji]. We use the 7-year WMAP data 
combined with large-scale structure, the Hubble constant 
measurement from the Hubble Space Telescope, super- 
novae type la, and big bang nucleosynthesis [1^ . We as- 
sume the flat A-Cold-Dark-Matter model with no fraction 
of massive neutrinos in the dark matter density (/^ = 0). 

In the likelihood analysis, we vary the following eight 
parameters: (i) baryon density today, fif,, (ii) dark mat- 
ter density today, Vic, (iii) the ratio of the sound horizon 
to the angular diameter distance, 9, (iv) the reionization 
optical depth, r;, (v) d{ko), (vi) ev-(fco), (vii) Vs {ko), and 
(viii) the Sunyaev-Zel'dovich amplitude, Agz- We take 




FIG. 2. 2D marginalized distribution for the quantum-gravity 
parameter 5{ko) = ao(5pi(A;o) and the slow-roll parameter 
eviko) with the pivot fco ~ 0.002 Mpc~^ for n — 2 and a — 2. 
The internal and external solid lines correspond to the 68% 
and 95% confidence levels, respectively. 



the pivot wave number feg = 0.002 Mpc~^ (£0 ~ 29) used 
by the WMAP team. 6{ko) and ev(fco) are constrained 
at this scale. While the bound on S depends on the pivot 
scale (and it tends to be smaller for larger ko), that on 
{koYS{ko) does not. 

While we assume a standard treatment of the reion- 
ization with a smooth interpolation, more general reion- 
ization scenarios can potentially affect constraints on ob- 



ll| . The analysis in 



servablcs especially for ?is > 1 
shows that the allowed region with rig < 1 is not strongly 
modified, which is the case for our potentials. 

The exponential term e~'^^ = {k^/kY in Eq. ([5]) gives 
rise to the enhancement of the power spectra on large 
scales, as we see in Fig. [T] In this sense the LQC cor- 
rections can be distinguished from the suppression effects 
coming from the noncommutative geometry or string cor- 
rections • For (T > 3, the growth of the term e~'^^ is 
so significant that 5p\{£) must be very much smaller than 
1 for most of the scales observed in the CMB, in order 
to satisfy the bound (5pi(£ = 2) ^ 1. More precisely, 
LQC corrections manifest themselves mainly at = 2, 3, 
where cosmic variance dominates, so it seems implausi- 
ble to isolate these effects. For cr < 3, the LQC modi- 
fication to the classical power spectra also affects larger 
multipoles t, and hence it is possible to constrain it from 
CMB anisotropies. 

In Fig. [2] we plot the 2D posterior distributions on the 
parameters 5{kQ) and ev(A:o) with fco = 0.002 Mpc~^ for 
n ~ 2 and cr = 2. The two parameters are constrained 
to be (5(fco) < 6.7 X lO'^ and tyiko) < 0.013 (95% C.L.). 
The modification of the large-scale power spectra (£ < 
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FIG. 3. 2D marginalized distributions as in Fig. (2] but for 
the case n — 2 and a = 1. 

20) shown in Fig. [1] leads to the upper bound on S{ko). 
The condition ^ gives the prior (5pi(£o) < 4.8 x 10"^ at 
£o = 29, so that for ao = 0(1) the observational bound 
is smaller by 2 orders of magnitude. For larger fcg the 
observational upper bounds on (S(fco) tend to be smaller 
for given a. For fcg = 0.05 Mpc~^ and ct = 2, we find 
that S{ko) < 1.2 x 10"'^ (95% C.L.), in which case the 
theoretical expected amplitude [Sp\{ko) ~ 10^^ or a few 
orders of magnitude higher 0]] can be accessible. 

For smaller a the observational upper bound on S{ko) 
tends to be larger, with milder enhancement of the power 
spectra on large scales. In Fig. [3] we show the likelihood 
results for cr = 1, in which case the LQC correction is 
constrained to be S{ko) < 3.6 x lO'^ (95% C.L.). Mean- 
while, the a-priori criterion ([6]) gives Sp\{ko) ^ 6.9x 10~^. 
For ao = 0{1), the case ct = 1 is marginally consistent 
with the combined SR/i5pi truncation. 

For cr < 1, the exponential factor e^'^^ does not change 
rapidly with smaller values of fs^t, so that the LQC effect 
on the power spectra would not be very significant even 
if S{kQ) was as large as ev{ko). Our likelihood analysis 
shows that the observational upper bound on 6{kQ) ex- 
ceeds the a-priori upper limit of 5pi(fco) given by Eq. ([6]). 
Since S{ko) can be as large as 1, the validity of the ap- 
proximation i5(A:o) < £5/(^0) used in the main formulas 
may break down in such cases. 

Under the conditions £ir,6 ^ 1, it follows that ~ 
(k^/2)((^///)^. Then the number of e-foldings during in- 
flation is given by iV = J^* di H fa k J^^ d(p j ^2(iv{^)-, 
where ipi is the field value at the end of inflation [de- 
termined by the condition ~ 0(1)]. For the power- 
law potentials one has N n/(4ev ) — n/4, which gives 
tv ~ n/{4N+n). For n = 2, the theoretically constrained 
range 45 < < 65 corresponds to 0.008 < < 0.011. 



The probability distributions of in Figs. [2] and [3] are 
consistent with this range even in the presence of the 
LQC corrections, so the quadratic potential is compati- 
ble with observations as in standard cosmology. 

In summary, in inflation combined with LQC inverse- 
volume corrections we provided general formulas for the 
scalar and tensor power spectra and placed observational 
bounds on the size of corrections for a quadratic poten- 
tial. In Q we ran the CosmoMC code also for other 
potentials such a.s V (x (p^ and V oc e"'*"'' (for which 
the inflationary observables reduce, again, to S and Cv). 
We found that the observational upper bounds are prac- 
tically independent of the infiaton potentials. This is 
because the LQC correction is approximately given by 
(5pi(fc) w 6pi{ko)(k/ko)~"' , which only depends on a and 
the pivot scale fep. Interesting and nontrivial effects do 
arise from the modified space-time structure underlying 
the dynamics. Even though quantum-geometry correc- 
tions are small, they can significantly change the run- 
nings of spectral indices. Thus, the observational bounds 
on (5pi can be much closer to theoretical values [O(10~^)] 
than often thought in quantum gravity. Our new tech- 
niques set the stage for systematic and stringent phe- 
nomcnological evaluations. 
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